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ABSTRACT 

The use of a self adaptive technique to compensate for the varying 
parameters of a flight control system is desirable for high performance 
aircraft and space vehicles. The technique discussed in this study con- 
Sists of using a sinusoidal signal to measure the change in the parameters 
and of compensating for the change by a gain adjustment, Analog and digital 
computer methods are used to investigate the performance of the system and 
the accuracy of the Amplitude Closed Loop Response Criterion for measuring 
damping ratio. The effects of varying airframe flexibility, center of 
gravity, frequency of the sinusoidal signal, and the Reference signal value 
are included. The applicability of the self adaptive technique for a 
particular vehicle depends onthe dynamic characteristics of the airframe 
and on the desired maximum variations of the performance characteristics. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Dr. G. J. Thaler, Professor of Electrical 
Engineering at the U. S. Naval Postgraduate School. The writer also 
wishes to express appreciation to Mr. L. K. Mattingly who supervised the 


study during a field trip at Autonetics in the summer of 1961, 
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CHAPTER I 


INTRODUCTION 


The dynamic characteristics of an airframe vary when changes in 
speed, altitude, position of center or gravity, and dynamic pressure 
occur, The variation in the dynamic characteristics cause variations 
in the response of the airframe to commands from a pilot or program un- 
less some method of compensation is used. Many autopilots employ some 
form of gain variation for compensation, For propeller and early jet 
aircraft the gain variation was accomplished by a schedule which was a 
function of some air data parameter such as dynamic pressure, altitude, 
or Mach number. In space vehicles and high-performance supersonic air- 
craft the gain scheduling becomes very complex or impossible for the 
following reasons: 

1. The dynamic and static characteristics of the airframe change 

rapidly over wide ranges. 

2, Air data measurements are not accurate or are not available at 

high speeds and altitudes, 

3. The airframe dynamic characteristics are not known with suf- 

ficient certainty to permit scheduling. 

To eliminate the requirement for scheduling and air data measurements, 
automatic control systems are being developed which include a self adap- 
tive loop. The purpose of the self adaptive loop is to maintain optimum 
response of the airframe by compensating for the changes in the dynamic 


characteristics without the necessity for external measurements. 








The Dither Seif-Adaptive System is one method for obtaining optimum 
response by self compensation, A study of this system for pitch control 
was undertaken with the following objectives: 

1. To determine the ability of the system to correct for the 

varying dynamic characteristics of a flexible airframe, 

2. To determine the optimum variation of the damping ratio and 
natural frequency of the short period mode during a typical 
flight and how this variation is affected by airframe bending 
and by changes in the center of gravity of the airframe, 

The first objective was investigated by the comparison of digital and 
analog computer results, The second objective was investigated primar- 
ily from an empirical approach although some analytical verification is 


included. 








CHAPTER 2 


DITHER SELF ADAPTIVE CONTROL SYSTEM 


2.1 Self Adaptive Systems 


Any closed loop control system has in effect some self-adaptive 
ability. In the usual control system a desirable static output is obtain- 
ed under varying input and open loop parameter conditions by proper design 
of the feedback loops, The name "self-adaptive' (1), however, is reserved 
for only those control systems which maintain optimum dynamic response 
under varying conditions instead of the static response that normal feed- 
back optimalizes. 

Self-adaptive systems can be classified according to the ability to 
automatically compensate for either changes in the system input or changes 
in the system parameters, such as environmental variations. The ideal self- 
adaptive system would compensate for both input and system parameter changes. 
Flight control systems are generally concerned with parameter changes and 
hence, are designed to compensate for only those changes. Since this in- 
vestigation is concerned with a flight control system, only the ability of the 
system to optimalize the dynamic response under varying aerodynamic para- 
meters will be considered. 

There are three basic operations which an adaptive loop must perform: 

1) A continuous measurement of system dynamic performance - The 
dynamic performance measuring method should have a negligible 
effect on the system response signal and, in like manner, the 
system response signal should have negligible effect on the 


measuring method. 








2) Á continuous evaluation of the dynamic performance on the 
basis of some predetermined criterion - The selection of the 
evaluation criterion depends on the physical capabilities of 
the system, Factors, such as realizability, complexity, and 
cost, may prevent incorporating into the system the ability to 
maintain a constant dynamic response, Compromises would have 
to be made so that the system would be physically capable of 
maintaining optimum a characteristic of the response; as, the 
relative stability - damping ratio - or the natural frequency. 
A figure of merit - which is a number or a method of compari- 
son - is then selected on the basis of accurately representing 
the selected dynamic response characteristic to be maintained 
constant, 

3) A continuous readjustment, based on the measured performance 
and evaluation, of system control parameters for optimum 
operation -- The methods used to adjust the parameters to 
obtain optimum operation can be listed under three headings: 
High Gain Linear Feedback, Programmed Compensation, and Compen- 
sation Using an Optimalizing Controller. 

In the High Gain Linear Feedback Method of adjusting performance under 
varying parameters the forward loop gain is maintained high by compensa- 
tion or with a relay. The desired dynamic resporse characteristics are 
contained in the feedback loop. The high gain with feedback makes the 
varying parameters in the forward loop negligible while the feedback loop 
determines the dynamic response of the system, This is not a new method and 


is used in many electronic applications for stabilization, 








he Programmed Lurpensatior Method measures the varying parameter, 

or the conditions which vary the parameter, and adjusts the compensating 
parameters in 4 controller in the forward loop to counteract the system 
variation, The adjustment is accomplished through a program which re- 
quires a knowledge of the relationship between the system and compensa- 
tion parameters. 

The Compensation Using an Optimalizing Controller Method is a type 
of feedback control which adjusts parameters so that a characteristic of 
the system response is optimalized; such as, the damping ratio. In the 
Controller Method an error signal using the output response level of the 
system can not be used as an error signal for the parameter adjustment as in 
a normal feedback method because the Controller is optimalizing a character- 
istic of the output and not necessarily the output level. The system may 
be used in conjunction with a model which gives the optimum characteristic 
in its output for a particular system input. The Controller samples the 
characteristic in the system output and compares the value with the optimum 
output of the model. The Controller then adjusts parameters until the model 


and system characteristic coincide. 


22 Description of Pitch Control System to be Investigated. 


The Dither Self-Adaptive System is capable of controlling a vehicle 
about its pitch, roll, and yaw axes. This investigation, however, will 
be concerned only with pitch control, The system which is to be investi- 
gated is shown in Figure 2,1. The position, rate and acceleration sensors 
that are required for the feedback loops are assumed to have a unity trans- 
fer function. The rectifier in the self adaptive loop is also assumed to 


have a unity transfer function. 
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In the conventional loop an error signal, d c, is formed by comparing 
the feedback signal to the command input signal. The error signal is 
amplified by the adaptive gain, ks , and the resulting signal operates the 
servo actuator. This investigation involves a missile which has gimbaled 
thrust chambers, The servo actuator controls the direction of the chamber 
axis which results in a thrust vector signal dt. The thrust vector signal 
affects the vehicle performance which is represented by the block labeled 
"Vehicle Dynamics". The reaction of the vehicle airframe to the change in 
direction of the thrust axis is detected by position, rate, and acceleration 
sensors and is fed back through the Inverse Model to the input. The function 
of the Inverse Model will be explained in Section 2.3. 

For the self adaptive loop a dither signal is introduced at the input. 
This signal is then amplified by kg , operates the servo actuator, and 
excites the vehicle dynamics. The pitch acceleration, 6, is passed through 
a bandpass filter where the dither frequency is separated from the command 
input signal frequencies. The dither frequency is then rectified and com- 
pared with a reference signal. The error, ESA) operates an integrating 
servo which varies the adaptive gain, kg و‎ until the reference and rectified 


Signals are equal. 


243 Principles Used to Maintain Optimum Dynamic Response. 


When a physical aerodynamic parameter changes, all the dynamic 
characteristic parameters of the system are affected. Typical physical 
aerodynamic parameters are Mach number and altitude, and typical dynamic 
characteristic parameters are natural frequencies and damping ratios. In 


terms of the system transfer function a change in aphysical parameter could 
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shif tke position ol all the poles and zeros of the transfer function. 
If constant dynamic respense was desired for the system, the variation in 
all the poles and zero would have to be counter-balanced. The complexity 
of the situation can be visualized by placing a root locus plot alongside 
a missile traveling many times the speed of sound. The physical parameters 
of the missile are constantly changing, and each parameter would be varying 
all the poles and zeros on the root locus. A typical root locus with vary- 
ing poles and zeros is shown in Figure 2.2. A controller would have to be 
built to sense each of the variations and then decide how to counter- 
balance by compensation. In practice exact compensation is too complicated 
to be useful and approximations are always made. 

Simplicity is one of the main objectives of the dither adaptive loop. 
The only components that have been added to the existing ptich control 
system are a dither generator, filter, rectifier, and a comparison-integra- 
tion network. No variable shaping networks are used, and all adjustable 
compensation is accomplished by varying the gain of the system. 

Aviators are primarily affected physically by the short period mode 
of the pitch aerodynamic equations. Tests have indicated that pilots prefer 
a damping ratio of 0.7 and a natural frequency of 3 radians per second for 
this mode (2). Since the missile used in this investigation is designed 
for manned flights, the damping ratio of 0.7 and frequency of 3 were de- 
sired. Since the dither system does not use variable shaping networks, it is 
physically impossible to maintain both the damping ratio and frequency con- 
stant at the desired values, Therefore, a damping ratio of 0.7 with a 
minimum natural frequency variation are defined as the optimum dynamic 


characteristics of the short period mode during a flight. 
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gure 2.2. Root Locus Pole and Zero Variation 
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There are actually two methods being employed in the overall system 
for varying the control parameter to obtain optimum performance. The 
control parameter of the system is the gain. The two methods, which are 
discussed in Section 2.1, are High Gain Linear Feedback and Programmed 
Compensation, 

If the gain of the system could be maintained at a high value, the 
desired dynamic response characteristic could be inserted in a feedback 
loop and would determine the dynamic response of the overall system. 
This is the principle of the High Gain Linear Feedback Method shown in 
Figure 2.3. .H(s) is referred to as the "Inverse Model" since the inverse 
of its transfer function determines the dynamic response of the system. 
If a damping ratio of 0.7 and a natural frequency of 3 radians per second 
were desired, the Inverse Model would have a transfer function of s o 
Unfortunately, the varying dynamic characteristics of the airframe, as 
represented by G2(s), cannot be submerged by merely using high gain feed- 
back loops (3). The required aerodynamic control moments and power required 
for high gains in flight control systems are not available due to limits 
imposed by vehicle design considerations unrelated to the automatic control 
System design. Even if high gains could be obtained, the gains could re- 
sult in severe bending and slosh mode oscillations when considering flexible 
airframes, The values that the Inverse Model transfer function can assume 
and still affect the short period mode in the desired manner are also limited 
for flexible airframes. The limit on the values of the transfer function 
for this investigation is discussed in Section 2.5. 

For flight control systems the limitation on high gains and on values 


of the Inverse Model prevents the use of a true High Gain Linear Feedback 
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self au-ptive scheme vori, the Tnverse Model can be used as a fixed 
shaping network to improve the dynamic characteristics of the short period 
mode and is employed as such in the system under investigation. 

In view of restricting the Inverse Model to the idea of a fixed shaping 
network, the self adaptive scheme used in the Dither System can be thought 
of as a form of the Programmed Compensation Method. The three basic opera- 
tions that an adaptive loop must perform are listed in Section 2.1. The 
manner in which the Dither System performs these operations are listed below: 

l. The dither signal is the means for the continuous measurement 
of the system dynamic response. In order to have a negligible 
effect on the system response signal, the dither frequency 
should be five or more times the short period mode frequency 
and the input level should be small so that the oscillations 
are not noticeable to the vehicle personnel. The reason for 
using pitch acceleration instead of pitch angle or rate as the 
pickoff for the adaptive loop is to enable smaller inputs for 
the same output. The upper limit for the dither frequency 
depends on the ability of the actuator to respond to the fre- 
quency and on the bandwidths of the rest of the system components. 
The upper limit is around 40 radians per second for available 
systems, 

2. The continuous evaluation of the dynamic performance is per- 
formed by comparing the amplitude of the pitch acceieration to 
dither input levei evaluated at the dither frequency to the 
Performance Criterion Reference, It remains to be discussed 
below and in Chapter 4 how accurately the dynamic performance 


is being evaluated, In this investigation the dampıng ratio 
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is being used as the measure of dynamic performance. 

3. The continuous readjustment of the system control parameter, 
which is the gain, is accomplished by a servo driven by the 
error resulting from the evaluation procedure, 

In the Programmed Compensation Method the parameter adjustment is ac- 
complished through a program which requires a knowledge of the relationship 
between the system and compensation parameters. In the dither system the 
program consists of a single reference level which requires a knowledge 
of the variation of the system damping ration with the amplitude of the 
acceleration to dither input level ratio. 

The degree of self adaptability of the dither system depends on the 
accuracy with which the damping ratio is measured by the acceleration to 
input level ratio. There are many figures of merit that can be used as a 
criterion for the dynamic performance of a system. Some of the figures of 

pO ጋ مه‎ መ pe 
merit are JE A Leat, 4 lel ZZ þet (Álit ner € 
0 
Ís an error signal. The applicability of these figures of merit for use 
in establishing a damping ratio of 0.7 has been determined for a second 
order system to be as follows: (4) 
مه‎ 
cU selects zero damping ratio as optimum, 
$e selects 0.5 damping ratio as optimum, 
00 
WP selects 0.7 damping ratio as optimum, 
Sea selects zero damping ratio as optimum, 
ነም. ማቸ 0.7 damping ratio as optimum. 
If the system was of second order, both gia and í el 
would select 0.7 as the optimum damping ratio based on the figure of 


merit approaching a minimum value at a damping ratio of 0.7. The 


oo 
ZE E core of merit has the sharpest minimum and would be the 
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preferred criterion, With a flexible airframe the system can not usually 
be approximated by a second order system with any accuracy, and the ap- 
plicability of the above criterion for establishing optimum damping ratics 
of the short period mode under these conditions is not known. The com- 

52 
ponents required to mechanize the D AE ER criterion into a 
system is certainly more complex than the Dither System, and the possibility 
of the criterion selecting a damping ratio of 0.7 for a non second order 
system is considered doubtful. A different figure of merit called the 
Impluse Response Area Ration (IRAR) will permit the selection of any de- 
sired damping ratio and have the error signal assume a zero value when the 
response is at the selected damping ratio (5). The method required for use 
of tbe IRAR criterion would be to ‘pulse the system and then measure the 
area of the response, The area is a measure of the damping ratio. By 
comparing the area with a reference area corresponding to the desired 
damping ratio an error signal is developed which would vary the gain of 
Ene Sistem until the error was zero. In a flexible airframe control system 
T does not use shaping networks, the poles and zeros of the vehicle . 
transfer function are continually varying. For the reference area to in- 
dicate the true damping ratio it would also have to vary with the parameters 
and in effect would need a self adaptive loop to adjust its value. The 
system would be very complex if not impossible to mechanize. 

The Dither System works on the principle of maintaining the gain of the 
system constant throughout the flight and assumes that a constant gain will 
maintain an optimum damping ratio. Referring to Figure 2.3 and substituting 

Em tine We AND OR = Op an adaptive loop is formed.‏ ےن 
T d^.‏ 
Sal) = Ks Ms GUF- Mo) CEK Uo).‏ 
Em (+ Ks M5 G, (346) 626) H(746)‏ 
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Since KsM, Gy (gus) Gi Yate) HG) has a denominator 
of greater order than the nunerator, £8 Gv.) = Ks Ms G, gx) CFG us) when 
Wo is a large value. Therefore, E (pu) is a direct measure of 
the gain, Ks Ms, of the system, The self adaptive loop maintains the 
value of Ks Mg equal to the reference value. If a constant value of 
KsMs is an accurate measure of the damping ratio, Y , for a flight 
control system, then a plot of Kg Mg versus £ should show a constant 
Ks Ms for a particular { for all flight cases. Using the Ms's of the 
flight cases to be investigated a plot of Ks Mg vs EM is shown in 
Figure 2.^ for aerodynamic equations with no bending modes included, the 
first bending mode included and the first and third bending modes included. 
The aerodynamic equations and conditions are discussed in Section 2.4. As 
is observed from the Figure, Ks Má is an accurate criterion for a con- 
stant 4 for no bending, but as bending modes are included and the order 
of the system increases KglVlg becomes a less accurate figure of merit. 
The investigation of the variation of the damping ratio for the various 
aerodynamic conditions and the determination of the optimum Performance 
Criterion Reference and the optimum dither frequency is contained in 


Chapter 4, 
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2.4 Aerodynamic Equations and Modifications. 

The airframe used in this investigation is the two staged Saturn 
Booster, Perturbated aerodynamic equations with coefficients for four 
flight cases are included in Appendix A. The flight cases are identified 
by the time in seconds after launch. The flight cases are t - 20, t - 75, 
t = 101.9, and t = 200. The t = 200 case is for the second stage only, 
as the first stage separates before t = 200 seconds. The reference for 
the equations is the nose of the second stage. For flight cases t = 20, t 
75, and t = 101.9 the acceleration, rate, and position sensors are located 
on the same platform in the first stage. The variation in performance for 
various platform positions will not be covered in this investigation. The 
platform position used had been selected for maximum system stability 
immediately after launch and for airframe considerations. The platform 
position for case t = 200 is in the nose of the second stage. The transfer 
functions which were determined from the equations using digital computers 
are included in Appendix A. 

The aerodynamic equations contain three bending modes and one fuel 
slosh mode. A typical root locus is shown in Figure 2.5, It will be 
observed from the Figure that the open loop system is unstable. Case t = 
20, which is the first case considered after launch and where the missile 
is still verticle, is the only flight case with a stable open loop. An 
unstable aerodynamic condition occurs when the center of pressure is for- 
ward of the center of gravity. A missile of conventional shape will be 
unstable unless stabilizing fins are attached because the center of 
pressure of the fuselage alone is close to the shoulder whereas the center 
of gravity is close to the geometric center (6). For large missiles, such 


as the Saturn, the trend is towards stabilization through a control system 
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rather than providing for a stable airframe because of the weight required 
for the stabilizing fins. The control torques required for the control 
system are provided for in the Saturn by gimbaled thrust chambers. It 
will also be observed from the Figure that the Slosh Mode has a negligible 
effect on the system stability. Multitank design with slosh suppressors 
are used in the Saturn to obtain this condition. The different slosh mode 
frequencies correspond to different tanks. 

The second bending mode becomes unstable for an adaptive variable gain 
of approximately 0.1 for the t = 20, 75, and 101.9 cases. Tn order to 
eliminate the possibility of obscuring the actual dither system performance 
by the complexity of adding bending mode cancellation schemes, the unstable 
second bending mode was handled in one of two ways: 

l. The second bending mode was eliminated from the aerodynamic 
equations assuming that a bending cancellation method was in use, 

2. The second bending mode was eliminated by assuming that a 
compensation method was in use, such as a frequency tracker and notch 
filter scheme, which placed a pole and zero over the open loop zero and 
pole respectively of the bending mode. 

A bending cancellation method which could be used for procedure 1 

O 
above is shown 6... ۴۶9۳ 2.6. This is one of many cancellation methods 
proposed by the Autonetics Division of North American Aviation. The purpose 
of the system is to form a signal of equal magnitude and opposite phase to 
the bending mode signal. The signal formed is then added to the control 
signal and effectively cancels the bending mode from the system response 
signal. Bandpass filters land 2 are variable tuned filters with center 


frequencies of Wo . If filter l is not tuned to the bending frequency, 
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the bending frequency -ignal component is phase shifted by an amount # 
and the output signal from filter 1 is sin (2/5 2 o The bending sig- 
nal is again phase shifted by an amount () by filter 2 with output signal 
being sin e+ a). The outputs from the two filters pass through 
phase shifters which produce a 90 degree phase shift between the signals. 
The shifted signals are then multiplied, and the resulting signal drives 
an integrator which varies the center frequencies of the filters until the 
output from the multiplier is zero. As is shown below, the output of the 
multiplier is zero when the filters are tuned to the bending frequency. 
Output of multiplier = sin 25200 cos Quo zd p) 

= (sin cos Y + cos WA? sin) (cos ህሬ.ሯ ር052ፖ -5፤። Wet sin 2D), 
The integrator is not sensitive to frequencies much greater than the 
bending frequency. 


Therefore, 


ዘ 


Effective Output 


1/2 cos 20 sin - 1/2 cos Y sin 27 
1/2 sing 


0 when Y is zero. 


With the filters tuned to the bending frequency, the 180 degree phase shift 
is accomplished by subtraction from the feedback signal. The amplitude of 
the bending signal is adjusted by the rectifier-integrator comparison loop 
using the outputs of filters 1 and 2, for the comparison. Another bend- 
ing mode cancellation method using a rejection filter philosophy is con- 
tained in Reference 7. 

The different aerodynamic conditions that are investigated are listed 
below: 

1. No bending and no slosh modes included in the equations. 


The resulting equations are the rigid body equations. 
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2. The tirst bending mode added to the rigid body equations. 

3. The first and third bending modes added to the rigid body 
equations, The second bending mode is assumed removed by a bending mode 
cancellation method. 

4. The first and third bending and first slosh modes added to 
the rigid body equations. The second bending mode is assumed removed by 
a frequency tracker and notch filter method,: 

In Chapter 3 the analog computer investigation was conducted using 
aerodynamic conditions 1, 2, and 3 and flight cases t - 20, 75, 101.9 and 
200. The second bending mode was included in condition 3 for flight case 
t = 200 because the mode is not unstable for this case. 

In Chapter 4, the empirical investigation was conducted using aero- 


dynamic conditions 1,2,3, and 4 and flight cases t = 20, 75, and 101.9. 


2.5 Limitations on Inverse Model by Flexible Airframe. 


The system under investigation is unstable unless a shaping network 
is employed. From a root locus viewpoint if a complex pair of zeros were 
provided in the left half plane, the unstable short period mode would be 
shaped into the left half plane and the system would be stabilized. The 
Inverse Mode! provides these shaping zeros. 

If a damping ratio of 0.7 and a natural frequency of 3 radian per 
second were the desired dynamic characteristics of the short period mode, 
the Inverse Model should ideally provide a pair of zeros at s = -2,1 د‎ 65 
on the root locus. A high system gain would then establish a complex pair 
of roots of the system characteristic equation at values s = -2.1+j 2.18. 
The roots would fix the dynamic characteristics of the short period mode 


at a damping ratio of 0.7 and a frequency of 3. As was discussed in Section 
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2,3, high garn capabilities ecr ust available in flight cortrol systems. 

The desired dynamic characteristics could still be obtained even 
though high gains are not available by increasing the real and imaginary 
coordinates of the Inverse Model zeros. The zeros could shape the short 
period dynamic characteristic curve in such a manner that for lower values 
of gain the desired characteristics are available. Increasing the coordin- 
ates of the Inverse Model zeros is possible in a rigid airframe with 
negligible center of gravity movement. For flexible airframes the bend- 
ing modes must be considered in establishing the position of the Inverse 
Model zeros. 

If the coordinates of the zeros are increased beyond certain values 
with a flexible airframe system, the zeros :hape primarily a different 
mode than the short period mode. While the short period mode would be 
stabilized, its maximum damping ratío and natural frequency values would 
be less and its coupling with the shaped mode would be greater than when 
the zeros shape primarily the short period mode. A digital computer investi- 
gation of the optimum position of the zeros was conducted. 

The preliminary investigation showed that flight case t = 75 was the 
critical flight case and that the two bending modes plus slosh mode condi- 
tion was the critical condition. Flight case t = 75 was further investi- 
gated as to the benefit of using rate and position feedback which would 
give one real axis zero or of using acceleration, rate, and position feed- 
back which would give the complex pair of zeros. The complex zero arrange- 
ment provided for more shaping of the short period mode and less possible 
coupling between the first bending mode and short period mode. The coordin- 


ates of the complex zeros were then varied to obtain optimum values. 
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Typical root loci are shown in Figure 2.7. For Inverse Model position 3, 
the zeros shape primarily the first bending mode. While the short period 
mode is stabilized, it is not being shaped in an optimum manner and has 
the disadvantages of decreasing in damping ratio as gain increases beyond 
a certain value and of being coupled with the bending mode. The short period 
mode could never obtain a 0.7 damping ratio. For Inverse Model positions 
l and 2, the short period mode is being primarily shaped. The bending mode 
is also shaped further into the left half plane increasing its stability 
while producing a negligible increase in coupling. Position 2, with the 
zeros at s = 1.8 +j.6, appears to be the optimum position. The maximum 
natural frequency of the short period mode at a damping ratio of 0.7 occurs 
at this position. As is observed, the natural frequency will never reach 
3 radians per second. A more complicated shaping network could be used to 
increase the natural frequency, but since one of the principal advantages 
of the Dither System is simplicity, a compromise has been made on the desired 
dynamic characteristics. As was defined in Section 2.3, the optimum dynamic 
characteristics of the short period mode for the investigation are a damp- 
ing ratio of 0.7 with a minimum natural frequency variation. 

The optimum Inverse Model transfer function for the severest flight 
case and condition will remain fixed at for all flight cases 


and for all aerodynamic conditions. 
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CHAPTER 3 


COMPUTER SIMULATION 


3.1 Introduction. 

The object of this Chapter is to investigate the ability of the 
Dither Self Adaptive System to correct for the varying dynamic character- 
istics of a flexible airframe during a typical flight. The correction is 
made by the adjustment of the control parameter which is the adaptive gain, 
kg 

An analog computer simulation was conducted for the primary purpose 
of determining the values of the adaptive gain which the system selected 
for the various flight cases and for the various aerodynamic conditions. 
The values of gain which were selected in the computer simulation are 
compared with the values of gain determined by digital computer methods. 
Digital computer techniques are employed in the investigation in Chapter 4, 
and the comparison with analog computer results in this Charter is used as 
an indication of the accuracy of the results of Charter 4 when the transi- 
tion is Eni from numerical to electrical simulation, The damping ratios 
of the short period modes which are associated with the gains selected by 
the analog simulation for the various aerodynamic conditions are determined 
from root loci using digital computer techniques and are compared to show 
the effect of a flexible airframe on the ability of the dither system to 
maintain a constant damping ratio. The damping ratios of the short period 
mode resulting from the use of the dither system are also compared with the 
damping ratios resulting from the use of a constant Kg system. 

The dither frequency and reference level which would maintain through 
the flight an average damping ratio of the short period mode of 0.7 with a 
minimum variation of the damping ratio and natural frequency were not used 
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in the cthputer sivulat on, Tre results of the computer sim.latien 
indicated the necessity ct determining optimum frequencies and reference 
levels. The investigation of this area is contained in Chapter 4. 

The secondary purpose of the analog computer simulation was to deter- 
mine for the various flight cases and aerodynamic conditions the rate at 
which the adaptive gain was adjusted and the reaction of the self adaptive 
loop to wind gusts, actuator movements, and loss of the dither signal. A 
complete analysis of the self adaptive loop was not conducted in the investi- 
gation; and therefore, the values of the loop gain and filter bandwidth for 
optimum dynamic response of the loop were not determined. The investigation 


was conducted using different loop gains for the various flight cases. A 


qualitative discussion of the adaptive loop is contained in Section 3.5. 


3.2 Analog Computer Setup. 


The analog investigation was conducted using aerodynamic conditions 
1,2, and 3 as described in Section 2.4, For Flight Case t = 200 the 
second bending mode was included because the mode is not unstable for 
this case. The system was simulated in two, forty amplifier, Pace Computers 
by Electronic Association Inc. The voltmeter used to determine the Kg 
values and the integrating servo and rectifier used in the adaptive loop 
were inherent components of the computers. Two six channel Sanborn re- 
corders were used. A Hewlett-Packard, Low Frequency, Function Generator was 
used as the source of the dither signal. 

The computer simulation diagram is shown in Figure B-1, The reference 
values of pitch acceleration, 8, pitch rate 0, and pitch angle, 8, which 
are measured at the nose of the second stage are converted to the platform 
values, Op , Op , and Op  , for flight cases t = 20, 75, and 101.9, be- 


cause the sensors for these cases are on the same platform which is located 


27 











in the first stage. Tne reference values are the platform values for 
Flight Case t = 200, because the platform for this case is in the nose of 
the second stage, 

The first flight case simulated was t = 75 for aerodynamic condition 
3. Static and dynamic checks were conducted, The procedure used for the 
dynamic checks was to excite the bending and short period modes by voltage 
steps introduced into the amplifiers of the various mode simulations and 
to compare the resulting frequencies and damping ratios with the frequencies 
and damping ratios determined by digital computer techniques. The trans- 
fer function resulting from the removal of the poles and zeros associated 
with the second bending mode from the vehicle dynamic transfer function, 
which was derived from the aerodynamic equations with three bending modes 
and one slosh mode, was considered a good approximation for determining the 
first and third bending mode frequencies for condition 3. A comparison with 
the analog results of the dynamic checks indicated that new transfer func- 
tions based on the applicable aerodynamic equations were required. The 
transfer function for flight case t = 75 for aerodynamic condition 3 in- 
dicated that the only poles and zeros of the transfer function, which was 
derived from the aerodynamic equations with three bending modes and one 
slosh mode, that were appreciably moved were the zeros associated with the 
third bending mode. The value of the zeros were changed from s --11.564j21. 
63 to s = -.17+ j24.3. The effect of the bending and slosh mode parameters 
on the vehicle dynamic characteristics can be determined by observing the 
movement of the poles and zeros of the transfer functions, which are contain- 
ed in Appendix A, for the various aerodynamic conditions. The static checks 


agreed to within one half of a per cent with analytical values, The mode 
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frequencies of the analog computer agreed to within two per cent of the 
frequencies determined by digital methods, and the damping ratios agreed 
to within five per cent. 

The procedure used to simulate the self adaptive loop was as follows: 

1. The filter was simulated and the bandwidth was varied until 
a sharp peaked frequency response was obtained with the peak at the dither 
frequency of 30 radians per second. A bandwidth of 3 radians per second 
produced a sharp peak, 

2. Two diodes which were inherent components of the computer 
were used to form a rectifier. The outputs from the two diodes were of 
different values; and therefore,gain pots were included with each of the 
diodes so that equal outputs could be obtained. 

3. The adaptive gain, kg , was simulated by using a linear 
gain pot driven by an integrating servo.’ The servo and pot were inherent 
components of the computer. One volt into the servo resulted in a gain 
change of 0.2, A linear pot was used because it was readily available. 

A shaped gain pot should be used in reality to obtain a small gain to servo 
input ratio for low values of gain and a large ratio for large values of 
gain. A small gain to servo input ratio is desired at low values of gain 
so that a wind gust or other disturbance would not reduce the gain apprec: - 
ably with the possibility of the system going unstable. As the gain in- 
creases the dynamics of the short period mode improve so that while a 
disturbance at high gains would vary the gain more than at low gains the 
response time is improved and stability considerations are not as critical, 
Compiomises must be made in the shaping of the gain pot between the large 
ratio of gain to servo input desired for fast adaptive response, the small 


ratio desired when disturbances are experienced, and the ratio required for 
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constant dynamic characteristics of the adaptive loop itself during 
a flight. The discussion of maintaining the adaptive loop gain appreciably 
constant during a flight is contained in Section 3.5. 

4. A gain pot was inserted before the integrating servo to 
provide for varying the ratio of Kg to servo input and the adaptive loop 
gain. The adaptive loop gain was adjusted by varying the gain pot, K:A ; 
until the gain adjustment response, Kg , ceased to oscillate and was 
critically damped with a suitable rate of change of the value of Kf 
The Reference DC signal was then adjusted to provide for a value of Kg 
of 1. The adjustment of the adaptive loop gain to give a suitable rate of 
change of Ks and the setting of the Reference level so that Kg had a value 
of 1 was performed for each aerodynamic condition of Flight Case t = 75 
in order to obtain reference values. The reference values that were obtain- 
‘ed eliminated the necessity for a complete evaluation of the analog simulat- 
ed op and were used in the other flight cases to form the basis for the 
comparison of the values of Kg that the system selected. 

When the system was energized, the servo chattered because it was 
following the half cycle sinusoidal like waveform at its input. A smooth- 
ing network would be required at the rectifier output to produce a flatter 
output and eliminate the servo chatter. In order to maintain the adaptive 
loop as simple as possible and also to investigate the effect of bending 
interference the bandwidth of the filter was increased to 5 radians. The 
third bending mode frequency which was now at the limits of the bandwidth 
was passing through the loop because thc waveform from the rectifier flat- 
tened out so that the servo stopped chattering. Since a dither frequency 
of 30 radians per second excited the third bending mode for aerodynamic 


condition 3, the computer results would indicate the performance of the 
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dith tem when bending or noise signals were contained in the filter- 
ed frequency of the self adaptive loop. In order to determine the per- 
formar» ur the system with no bending and minimum noise signals in the 
self adaptive loop, reference values were also recorded for Flight Case 
t = 75 for aerodynamic condition 3 with a dither frequency of 20 radians 
per second. The filter was adjusted for a center frequency of 20 and a 
bandwidth of 3. The lower frequency did not cause the servo to chatter. 
The general procedure used in the simulation for Cases t = 20, 101.9, 

and 200 was as follows: 

1, The computer was set up for condition 3 of the particular 
flight case. 

2, Static and dynamic checks were conducted. 

3. The Reference level and KsA were set at the values used 


for the corresponding aerodynamic condition and dither frequency of Flight 


^, The values of Kg as read from the computer voltmeter were 
recorded for aerodynamic conditions 1, 2, and 3 with a dither frequency 
of 30 radians per second and for condition 3 with a dither frequency of 
20 radians per second, 

5. The value of KsA was adjusted to give a suitable rate of 
change of Kg for the particular aerodynamic condition and flight case. 
The Reference level was then adjusted to give the corresponding value of 
Kg as determined in Step 4 with a dither frequency of 30. 

6. Recordings were made at the dither frequency of 30 for 
aerodynamic conditions 1 and 3. During the run for each condition, Kg 
was offset by a step voltage put into the integrating servo amplifier; an 


actuator signal, gy , Was put into the system by a step voltage into the 
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simulation for the actuator; a wind gust was simulated by a step voltage 
into the attack angle, Of , amplifier; and a dither signal failure was 
simulated by turning off the signal generator. 


Step 6 was also conducted for Case t = 75. 


3.3 Comparison of System Selected Gains. 


The values of the adaptive gain, Kg , which were recorded by the 
procedures described in Section 3.2 during the analog computer simulation 
are compared in this Section with values of gain determined from digital 
methods to establish a correspondnece between the two procedures, A Recomp 
IL Computer was used for the digital calculations. The damping ratios of 
the short period modes corresponding to the recorded adaptive gains are 
compared to establish the relative adaptibility of the system under vary- 
ing aerodynamic conditions and dither frequencies and are compared with 
the damping ratios corresponding to a fixed gain system to establish the 
benefit of the adaptive system for the airframe being used for the investi- 
gation, The analog computer simulation was conducted for aerodynamic 
conditions 1, Z, and 3 as described in Section 2.4. The digital investiga- 
tion included aerodynamic condition 4, and the damping ratios associated 
with the gains for a particular Dither Reference level for this aerodynamic 
condition are included in the comparison of the damping ratios of the ana- 
log computer. 

The determination by digital methods of the values of Kg which the 
system would select was accomplished by computing the amplitudes of the 
closed loop response for various gains and for frequencies of 20 and 30 
radians per second. Plots of the amplitudes of the closed loop response 


versus Ks for a dither frequency of 30 were made of the flight cases for 








each aerodynamic condition. For aerodynamic condition 3 plots were made 
for dither frequencies of 30 and 20. The plots are shown on Figures B-2 
(a) through B-2(f). An adaptive gain of 1 was set on the analog computer 
for each aerodynamic condition of Flight Case t = 75 and was used as the 
standard for determining the values of Reference level and adaptive loop 
gain pot settings which were maintained constant in recording the adaptive 
gains of the remaining flight cases, The adaptive gain value of 1 for Case 
t = 75 was also used as the standard for the digital method. The reference 
amplitude of the system was determined on the plots by the intersection of 
the Kg equal to 1 line with the Flight Case t = 75 curve. The reference 
amplitude eliminated the necessity of evaluating the components of the 
adaptive loop to determine a Reference Signal level and was used to deter- 
mine the values of gain which the system would select for the other flight 
cases by the intersection of the reference amplitude and the particular 
flight case curve. 

The comparison of the adaptive gains determined by analog and digital 
methods is shown in Table I and II for aerodynamic conditions land 2 
respectively at a dither frequency of 30. The reference amplitude deter 
mined for aerodynamic condition 3 at a dither frequency of 30 was too large 
for use with Flight Cases t - 101.9 & 200 as shown on Figures B-2(c). The 
dither frequency was varied in the digital method for each flight case 
until a reference amplitude was determined which resulted in the gains which 
were selected by the analog computer. The dither frequencies required to 
obtain the gains are shown in Table III. The plots of closed loop response 
amplitude versus Kg for the different dither frequencies required are shown 
in Figure B-2(d). The third bending mode was excited at a dither frequency 


of 30 radians per second and the filter bandwidth allowed the bending signal 
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| TABLE I 


Gain Comparison for Aerodynamic Condition 1 at Dither Frequency of 30 








| Flight Case | Gain Selected by | Gain Determined by | 
Analog Computer Digital Procedure 
Bu 20 1 | 1.69 چ‎ 
t = 75 نی‎ | 1.) 
10 101,9 0.54 07333 
2:200 IE 5 +17 


TABLE II 


Gain Comparison for Aerodynamic Condition 2 at Dither Frequency of 30 





ت سو RETTET‏ 


Flight Case | Gain Selected by Gain Determined by 
Analog Computer Digital Procedure 
pn 20 1 1.3 
= 101,9 0.8 0.72 


Aoro ې‎ 


t = 200 3.4 3 7 | 
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to pass thrcugh the adaptive luop as was discussed in Section 3.2. 

The dither signal source could not be accurately set to within 0.628 

radians per second; and since two dither frequencies were used for each 
flight case, the signal source had to be adjusted for each case. The 
variation in the digital frequencies required to obtain the analog gains 

is approximately within the accuracy of the signal source settings. The 
necessity for very accurate determination and setting of the dither fre- 
quency and Reference Signal level when bending modes are excited close to 

the dither frequency is illustrated by the difficulty in obtaining a compari- 
son of the gains for the digital and analog methods. When bending or noise 
signals pass through the adaptive loop, there is a possibility of the system 
having a choice of selecting either of two adaptive gains which are relative- 
ly close in magnitude. The selection of the wrong gain might cause the 
system to be unstable. As is shown in Figure B-2(d), the reference ampli- 
tude intersects the flight cases at two relatively close values of Kg . 

The comparison of the gains determined for aerodynamic condition 3 at a 
dither frequency of 20 is shown in Table IV. The difficulties experienced 
with a dither frequency of 30 were not experienced at a frequency of 20 
where bending modes are not excited. The plots of closed loop response 
versus Ke for a dither frequency of 20 are shown in Figure B-2(e). 

Root loci were determined for the various flight cases and aerodynamic 
conditions using a Recomp II digital computer. The damping ratios of the 
short period modes corresponding to the gains selected by the analog computer 
with a dither frequency of 30 and to the gains selected by digital techniques 
for aerodynamic condition 4 with a dither frequency of 30 were determined 


from the root loci and are compared in Table V. 
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TABLE III 


Frequencies Required by Digital Method for Grins Determined by Analog 
Method for Aerodynamic Condition 3 at Dither Frequency of 30 






Flight Case | Gain Selected by 
Analog Computer 





Frequency Required 
by Digital Method 
in Radians per Second 





| t = 20 1.0 

$ 

mee = 75 09 28.5 
EN 9 ۱ 2] 30,0 
1 : 


t = 200 5.0 28.17 | 


TABLE IV 


Gain Comparison for Aerodynamic Condition 3 at Dither Frequency of 20 
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TABLETY 


Damping Ratio Comparison for Flights Under Various Aerodynamic Conditions 
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Flight Case t دم‎ was used as the reference for both the analog and 
digital procedures; and theretore, the maximum variations of the damping 
ratios of a particular flight from the damping ratio of cases t = 75 is 
probably the best criterion co show the effect of bending modes on the 
adaptibility of the system, The results of the comparison of the first 
three aerodynamic conditions shows adaptibility decreasing as bending is 
included, but the comparison is valid only for a dither frequency of 30 
and must be checked for other dither frequencies. The investigation as 
to the effect on the results of a comparison when the dither frequency is 
varied is contained in Chapter 4. A different method was used in aero- 
dynamic condition 4 than in condition 3 to remove the second bending mode. 
The different methods resulted in a large variation of the position of 
the zeros associated with the third bending mode as is noted by examining 
the transfer functions in Appendix A. The comparison of the two conditions 
in Table V indicates that for a dither frequency of 30 the frequency track- 
ing and notch filter scheme for removing the second bending mode results in 
a more adaptive system than the bending cancellation scheme. The total varia- 
tion of the damping ratios of the two conditions is the same but condition 4 
includes a slosh mode, has the larger value of variation from the reference 
value of the t = 75 case towards greater damping, and has an average damping 
ratio during the flight less than condition 3. For a larger average damping 
ratio during a flight, a smaller variation in the damping ratios would be 
expected because a larger average damping ratio means larger gains are used 
and large gains force the damping ratios towards the constant damping ratio 
determined by the position of the Inverse Model. 

The damping ratios as determined from the root loci for the gains select- 


ed by the analog computer for aerodynamic condition 3 with dither frequencies 


38 








of 30 ati 20 dec € im rable Vi with the damping ratios of a fixed | 
gain system vith the J410 sel at a value of 1. The gain of a value of | 

was selected for the fixed gain system in order to establish the same 

reference for Case t = 75 that was used for the Dither System. The results 

of the comparison show that the fixed gain system is more adaptive than 

the dither system at a frequency of 30 radians per second but is less 

adaptive than the dither system at a frequency of 20. The great difference 

in the performance of the dither system as the dither frequency was varied 

led to the investigation of obtaining the optimum dither frequencies for 

the different aerodynamic conditions. This investigation is contained in 


Chapter 4. 


3.4 Reaction of the System to Disturbances. 


Aerodynamic conditions 1 and 3 were used in the investigation of the 
rate at which the adaptive gain was adjusted and the reaction of the adap- 
tive loop to wind gusts, actuator movements, and loss of dither signal. 

The investigation was conducted with a dither frequency of 30 radians per 
second so that the evaulation of the adaptive loop for aerodynamic condition 
3 would be performed under the unfavorable conditions of having bending 
signals included with the filtered dither signal in the adaptive loop. 
Since a complete analysis of the adaptive loop was not performed, the 
adaptive loop gain was adjusted for each flight case with the following 
considerations: 

l. To produce a suitable rate of change of the adaptive gain, Kg : 

2. To show the variation in the loop performance for different 
flight cases with the same setting of the servo gain pot, KgA . 

3. To use several values of ریک‎ for the flight cases with one 


value producing an overshoot in the Ks response. 
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Damping Ratio Comparison of Dither System and Fixed Gain System 


est^  .መው 





tempor tal? s, I 














TABLE VI 
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The 51 1 ۲ abed for corresponding flight cases under the 
Jiffexe TE gui. COplittons to show the effect of bending on the 
loop perforuoti, 

16601:1111 3 were wade of the reaction of the system to various distur- 
bances for the different flight cases and aerodynami. conditions, The traces 
for the flight cases for aerodynamic condition 3 showing the actuator 
command signal dy; the attack angle, & ; the adaptive gain response, Kg; 
the signai from the adaptive loop rectifier, SAE; and the first bending 


mode, d,, are contained in Figures B-3(a) through B-3(d). The procedure 


1? 
used in obtaining the traces is contained in Steps 5 and 6 of the procedure 
used to simulate Cases t = 20, 101.9, and 200 in Section 3.2. 

The comparison of the rates of change of Kg when Kg is offset above 
and below the undisturbed magnitudes selected by the system for the vari- 
ous cases and conditions is shown in Table VII. Cases t = 20 and 75 have 
the same Keg setting, or Kg to servo input ratio, but th. response times 
vary considerably. The variaticn illustrates the fact that the Kg to 
servo input ratio must also be programmed for the flight if constant dynamic 
conditicns are to exist in the system. There must be an optimum Reference 
Level, Dither Frequency, and Adaptive Loop Gain for optimum operation of the 
dither system. The program required to maintain the optimum value of the 
Kg to servo input ratio can be placed into the loop in the form of a shaped 
gain curve for pot هی‎ as is explained in Section 3.5. The adaptive loop 
gain for Case t = 101.9 resulted in an overshoot in the Kg response as is 
shown on Figure B-3(c) The response time for this case is the fastest as 
would be expected since the Kg responses for the other cases are ov. damped. 
As bending was added the response time increased except for Case t = 101.9 


where the response time remained constant. In order to form some basis 
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Comparison vi Kates of Change of the Adaptive Gain for Various Loop Gains 
and Aerodynamic Conditions 
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for establishing ہ11‎ the times listed in Table VII are sufficiently fast for 
the system to function properly, the time intervals between flight cases 
are divided by the gain changes required between flight cases to estab- 
lish figures of merit. The figtres of merit are based on the assumption 
that Kg is changing linearly with the time between flight cases and are 
listed below: 

Eu 75 OOseconds per unitychange in Kg 

== ያን. ሽን) 24,6 seconds per unity change in Kg 

tx C101:9 0507200 33.8 seconds per unity change in Kg 
Even though the loop was not optimized the times listed in Table VII are 
well within the values established as figures of merit. 

The comparison of the reaction of the adaptive loop to actuator signals 
for the various flight cases and aerodynamic conditions is shown in Table 
VIII. The magnitude of Kg was changed in all cases and conditions when an 
actuator signal was introduced into the system, but the adaptive loop re- 
turned ا‎ , while the signal was still being applied, to the original 
magnitudes for all cases and conditions in relatively small time intervals. 
The ideal system would result in no change in KS when an actuator signal 
was experienced. Excessive per cent changes in Kg occurred for Flight 
Case t = 20 for aerodynamic conditions 1 and 3 and for Flight Case t = 
101.9 for aerodynamic condition 3, If the adaptive loop had an optimized 
gain with a properly shaped Kg gain pot instead of the linear pot used in 
the simulation, these large changes in Kg would be greatly reduced. One 
of the considerations that must be used in arriving at the proper shape for 
the gain curve of the pot is to have a small ratio of Kg to servo input for 
small values of Kg, The undisturbed magnitudes of Kg for Case t = 20 was 


l and for t = 101.9 was 2.1 so that the shaped pot would aid these cases 
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Comparison of the Reaction of the Adaptive Loop to 
Actuator Signals 
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in preventirg lacge chádeco in ke when disturbances were experienced, 
The same magnitude actuator signals for the bending aerodynamic condition 
produced approximately the same change in the angle of attack as for the 
no bending aerodynamic condition, but the bending condition resulted in 
greater changes in Kg and longer time intervals to return to conditions 
that existed before the disturbance than the no bending condition. 

The comparison of the reaction of the adaptive loop to wind gusts for 
the various flight cases and aerodynamic conditions is shown in Table IX. 
A one degree change in the angle of attack corresponds to wind velocities 
greater than would be expected by the missile during the flight. The per 
cent change in Ks for all flight cases and conditions is small, within 1 
per cent. The values of Ks are returned to the undisturbed magnitudes 
while the gust is still being experienced for all flight cases and condi- 
tions in small intervals of time. The per cent change in Ks and the 
time interval to return to conditions that existed before the disturbance 
increase as the airframe becomes more flexible, 

The Reference signal is set so that if the dither signal source fails 
the value of Kg is increased to insure that the system will not become 
unstable. The comparison of the rates of change of Kg for the various 
flight cases and aerodynamic conditions when a dither failure occurs is 
shown in Table X. The rate of increase of Kr decreases as the airframe 
becomes more flexible. 

The results of the investigation contained in this Chapter show that 
for the computer simulation of the dither system the performance of the 
adaptive loop deteriorated when the airframe became flexible, The adaptive 
gain is varied by actuator signals and wind gusts while the ideal system 


would have no variation of gain when experiencing disturbances. The 
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variations of gain for wind gusts were small and were corrected during 
the gusts ina relatively short period of time. The variations of gain 
for actuator signals were large for two flight cases and small for the 
remaining flight cases. The variations of gain were corrected during 
the signal in relatively short periods of time, The rates of change of 
adaptive gain were satisfacto: ly based on the criterion of a constant 
rate of change of gain required between flight cases. By optimizing 
the gain of the self adaptive loop and properly shaping the adaptive 
gain pot the rates of change of gain could be increased, the variation 
in the adaptive gain due to actuator and wind gusts could be greatly 
reduced, and the time required to correct for disturbances could be 


reduced. 


3.5 Self Adaptive Loop. 


The importance of optimizing the gain of the adaptive loop was shown 
in Section 3.4. For small perturbations the loop is linearized in order 
that a qualitative investigation may be conducted (8). The linearized 
loop with the Reference value as the input is shown in Figure 3.1. The 
Reference value is compared with the rectified filtered dither signal, 
and an error signal is formed, The error signal crives the integrating 
servo with the output shaft positic- indicated by 2 . The adaptive gain 
pot has a transfer function of ha which is the slope of the curve of 
Ks versus servo shaft position. The gain, Ks و‎ is then multiplied by 

.. 
the vehicle dynamics transfer function, a Ks و‎ which is the slope of 
the Closed Loop Response Amplitude versus Kg curve at the Reference 
value for the particular flight under consideration. The resulting signal, 


ود 


0, ís then fed back through the filter. The characteristic equation of the 
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loop "ASE = - l. The root locus is shewn in Figure 3.1. 

As the loop gain increases, the response oscillates. It was found dur- 
ing the analog computer study that if the gain of the adaptive loop was 
increased beyond a certain value for a particular flight case that the 

Ks response would oscillate and as the gain was reduced the Kg response 
became over damped. The observations of the study agree with the root 
locus. 

In order to maintain the adaptive loop gain constant at the optimum 
value that is determined, the product of the Adaptive Gain Pot transfer 
function and the Vehicle Dynamic transfer function must be a constant. If 
the pot is shaped so that 5ታ is the reciprocal of 7 at the value 
of Ke that the system selects for the particular flight case, the value 
of Kv will be a constant equal to 1. The shape of the gain curve for the 
pot for aerodynamic condition 3 with a dither frequency of 30 radians per 
second is shown in Figure 3.2. 

In an actual design a fixed loop gain could be established by varying 
KA and Kg until a desirable time response was obtained for the servo out- 
put shaft, and the gain of the Kg pot could be shaped as explained above. 
If the performance of the loop was not satisfactory as to the rate of gain 
change, the change in gain caused by disturbances, or the time required to 
correct the gain after a disturbance, compromises would have to be made. 
The slope of the gain curve could be veried to improve the performance. 

In varying the slope, however, Kv would vary and the adaptive loop gain would 
no longer be a constant for a flight. A varying Kv would vary the time 
response of the servo shaft, and the performance of the loop would be af- 
fected. A different dither frequency with its corresponding Reference 

value could vary the shape of the gain curve resulting in the satisfactory 
performance of the loop. Optimum setting of the loop gains and shape of the 
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CHAPTER á 


CLOSED LOOP RESPONSE AMPLITUDE AS A CRITERION FOR CONSTANT DYNAMIC RESPONSE 


4,1 Introduction. 

The object of this Chapter is to determine for the dither system the 
accuracy of evaluating the dynamic performance of a variable parameter con- 
trol system - a flight control system - by the amplitude of the closed loop 
response, A dither signal is used as an input to the system and is varied 
by the dynamic characteristics of the system. The dynamic characteristics 
are then evaluated by the comparison of the amplitude of the dither signal 
with a fixed level Reference signal which represents the desired dynamic 
characteristics, The optimum dynamic performance of the system under in- 
vestigation was defined in Section 2.3 as a damping ratio of the short 
period mode of 0.7 with a minimum variation of the natural frequency. 
Therefore, the accuracy of evaluating the dynamic performance by the ampli- 
tude of the closed loop response is based primarily on the accuracy of 
evaluating the damping ratio of 0.7. 

The Reference signal values for dither frequencies of 20 and 30 
radians per second were not determined for an average damping ratio of 0.7 
for the flight for the computer sim lztion discussed in Chapter 3. The 
Reference values used, however, resulted in approximately equal average 
damping ratios for the flights, but the maximum variation of the damping 
ratios was twice as large for the dither frequency of 30 as for the frequency 
of 20. The maximum variation of the damping ratios of a flight also differed 
with the flexibility of the airframe and with the consideration of variable 
fuel load. The effect of varying the dither frequency, of increasing the 
flexibility of the airframe, and of varying the center of gravity on the 
accuracy of the performance evaluation will be determined in this Chapter. 


The results of the investigation will permit the selection of the Reference 
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signal value and the dither frequency, or frequencies, which result in 
the optimum variation of the damping ratio and natural frequency of the 
short period mode from the desired system dynamic performance. 

The investigation was conducted principally from an empirical view- 
point although some analytical verification is included in Appendix C. 
Aerodynamic Conditions 1, 2, 3, and 4 and Flight Cases t = 20, 75, and 
101.9, which are described in Section 2.4, are included in the investigation. 
Flight Case t = 200 was not included because the sensors are located ina 
different location than the sensors for the other cases. The performance of 
flights evaluated from cases that are measuring the reactions of the system 
at different positions in the airframe was not considered a true indication 
of the capabilities of the dither system. Therefore, the performance of 
flights are evaluated in this investigation with only flight cases that 
are measuring reactions at the same airframe position. Application of 
the results to a general variable parameter system is discussed in Section 


a4. 


4.2 Procedure, 

The investigation of the effect of varying the dither frequency on the 
accuracy of the closed loop response amplitude in determining the damping 
ratio of 0.7 was conducted using Cases t = 20, 75, and 101.9 as representa- 
tive of a typical flight. The dither frequency was varied in the range 
from O to 40 radians per second. The upper limit of 40 radians was select- 
ed because available actuators will not respond to higher frequencies. For 
various frequencies in the dither range of each of the four aerodynamic 
conditions Reference signal values were determined which resulted in an 
average damping ratio during the flight of 0.7. At the Reference value for 
a particular aerodynamic condition and dither frequency the maximum varia- 
tion of the damping ratios during the flight was determined. Plots of the 
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maximum damping ratio variation, AS , versus the dither frequency, W , 
were then made for the four aerodynamic conditions. The procedure employ- 
ed in determining the Reference value and the maximum damping ratio varia- 
tion will be described by the use of a typical example. A dither frequency 
of 22 radians per second and aerodynamic condition 3 will be used for the 
example. 

Root loci for the three flight cases were determined for aerodynamic 
condition 3 using a Recomp II digital computer. The damping ratios, کے‎ ; 
of the short period modes associated with the adaptive gains, Kf , were 
determined from the root loci, and plots of § versus Kg were made for the 
flight cases as shown in Figure 4.1. The amplitudes of the condition 3 
closed loop transfer functions were evaluated at the dither frequency of 
22 radians per second for various values of adaptive gain. The plots of 
closed loop response amplitude versus adaptive gain for the flight cases 
are shown in Figure 4.2, A graph of the damping ratio versus flight case 
was constructed as shown in Figure 4.3. The procedure used in establishing 
the closed loop response amplitude which would result in an average damping 
ratio during the flight of 0.7 was as follows: 

1. An amplitude for the closed loop response was selected, and 
the values of K$ corresponding to the amplitude were determined for the 
three flight cases from the Amplitude versus Kf plots. 

2. The values of the damping ratios corresponding to the values 
of Ka, which resulted from the selected amplitude, were determined from 
the d versus Kg plots. 

3. The values of damping ratios were plotted on the damping ratio 
versus flight case graph. 

4, If the average of the three damping ratios corresponding to 
the same closed loop response amplitude was not 0.7, a different response 


amplitude was selected. Steps 1 through 4 were repeated until the average 
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damping ratio was 0.7. 

When the closed loop response amplitude which resulted in an average 
damping ratio of 0.7 was determined, the difference between the maximum 
and minimum values of the damping ratios for the three flight cases was 
taken as the maximum variation of the damping ratio, A$ , during the 
flight for a dither frequency of 22 radians per second. The maximum 
variation and dither frequency were plotted on the AS versus W graph, 
shown in Figure 4.4, for aerodynamic condition 3. The procedure was 
repeated for various frequencies for the various aerodynamic conditions 
to construct the AS versus Wplots for the dither frequency range. 

The optimum dynamic performance was defined as a damping ratio of 0,7 
with a minimum variation of the natural frequency. The procedure that has 
been described determined the variation of the damping ratios with varying 
dither frequencies while maintaining an average ratio of 0,7 during a flight. 
The procedure that will be described determined the variation of the natural 
frequencies of the short period modes with varying dither frequencies while 
maintaining the same average damping ratio as that used in determining the 
damping ratio variation. The reference for determining the damping ratio 
variation was the defined optimum performance ratio of 0.7. The reference 
for determining the natural frequency variation was selected as the average 
natural frequency of the short period modes of the three flight cases at the 
damping ratio of 0.7. The procedure employed in determining the natural 
frequency variation will be described by the use of the same example that 
was used for the damping ratio variation: namely, a dither frequency of 22 
radians per second and aer: dynamic condition 3. 

The natural frequencies of the short period modes at a damping ratio 


of 0.7 were determined from the root loci for the three flight cases. 
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The three frequencies were averaged to determine the reference natural 
frequency for condition 3. The associated adaptive gains, Kg , and 
natural frequencies, Wn, of the short period modes were determined from 
the root loci, and plots of Wo versus Kg were made for the flight cases 
as are shown in Figure 4.5. The procedure used to determine the frequency 
variation was as follows: 

1, The closed loop response amplitude which resulted in an 
average damping ratio of 0.7 for the flight was determined from the damping 
ratio versus flight case graph for the dither frequency of 22, 

2. The adaptive gains for the flight cases which correspond to 
the amplitude which resulted in the average damping ratio of 0.7 were deter- 
mined from the Amplitude versus Kç plots. 

3. The adaptive gains determined from the Amplitude versus Kp 
plots were used to enter the W versus Kf plots to determine the natural 
frequencies of the three flight cases. 

4. The three natural frequencies were averaged, The reference 
natural frequency was then subtracted from the averaged frequency to 
determine the natural frequency variation for the dither frequency of 22 
radians per second. 

The natural frequency variation, JIT was plotted on the AU 
versus 2L/ graph which is shown in Figure 4-6 for aerodynamic condition 3. 
The procedure was repeated for various frequencies of the various aerodynamic 
conditions to construct the AW versus Wplots for the dither frequency range. 

The plots of AÇ versus W and AW, versus W for the four aerodynamic 
conditions will be used in Section 4.3 to determine the accuracy of the 
closed loop response amplitude as a criterion for constant dynamic response 
and the effect on the accuracy of varying the dither frequencies, increasing 


the flexibility of the airframe, and varying the center of gravity. 
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4.3 Results of the Empirical Investigation. 





One of the basic operations, which are discuss. d in Section 2,1, 
of an adaptive loop is the continuous measurement of che system dynamic 
performance. The measuring method must have a negligible effect on the 
system response signal. In the dither self adaptive flight control system 
the measuring metthod involves a dither signal, and the important system 
response signal is the short period mode response. The possible frequencies 
of the short period mode response are never greater than 1.5 radians per 
second for the cases and aerodynamic conditions considered. To insure 
that the dither signal does not interfere with the short period mode re- 
sponse the lower limit for the dither frequency range was set at 10 radians 
per second. The upper limit was previously set at 40 radians per second 
to conform to the characteristics of available actuators. 

The effect of varying the dither frequency, increasing the flexibility 
of the airframe, and varying the center of gravity on the accuracy of evalu- 
ating the specified damping ratio is shown in Figure 4.4, The AS versus 

W plots are based on the requirement that an average damping ratio of 0.7 
is maintained for the flight. The accuracy with which the damping ratio is 
evaluated for aerodynamic conditions 1 and 2 is independent of the dither 
frequency in the range from 10 to 40 radians per second. Ihe damping ratio 
for aerodynamic condition 1, which is the rigid airframe condition, is very 
accurately evaluated with a variation of 0.06 between the high and low values 
of the damping ratios during the flight. As the airframe becomes flexible 
by including the first bending mode, aerodynamic condition 2, the accuracy 
decreases, The variation of the damping ratios for condition 2 is approxi- 
mately 0.15. As the flexibility of the airframe increases by including bend- 


ing modes 1, 2, and 3 but eliminating mode 2 by a bending mode cancellation 
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method, aerodynamic condition 3, the accuracy of evaluation the damping 
ratio depends on the dither frequency. A difference in the variation of 
the damping ratios for dither frequencies of 20 and 30 radians per second 
was determined in the computer simulation, discussed in Chapter 3, for 
this aerodynamic condition. The average variation of the damping ratios 
for the dither frequency range is approximately 0.24 with a minimum varia- 
tion at a dither frequency of 23 radians per second of approximately the 
same magnitude as that for a rigid airframe: namely, 0.06. If bending 
modes 1, 2, and 3 with bending mode 2 eliminated by a frequency tracking 
and notch filter method and a varying center of gravity in the form of a 
slosh mode are included in the airframe equations, the accuracy of evalua- 
ting the damping ratio is independent of the dither frequency from 12 to 
35 radians per second. The variation of the damping ratios for aerodynamic 
condition 4 in the 12 to 35 frequency range is approximately 0.435. For 
dither frequencies greater than 35 the variation decreases and is 0.375 at 
a frequency of 40 radians per second. The different methods used in eliminat- 
ing the second bending mode in aerodynamic conditions 3 and 4 resulted in 
large differences in the values of the real parts of the zeros associated 
with the third bending mode of the airframe transfer function as is discussed 
in Section 3.2. The zero for the third bending mode for condition 3 is close 
to the imaginary axis while the zero for condition á is not. 

In summary, the frequency in the dither range has no appreciable affect 
on the accuracy of evaluating the damping ratio for aerodynamic conditions 
ہلا‎ 2, and 4. The accuracy depends greatly on the choice of the dither fre- 
quency for condition 3 with the magnitude of the variation varying from the 
small values obtained in condition 1 to the large values obtained in condi- 


tion 4. The accuracy decreases as the airframe becomes flexible and as a 
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varying center of gravity is included, The large variations of condi- 
tion 4 could be reduced Fy using an average damping ratio of greater 

than 0.7 for the flight because for large damping ratios the amount of 
change of gain per unit change in damping ratio is greater than for small 
damping ratios as is shown in Figure 4.1. 

The effect of varying the dither frequency, increasing the flexibility 
of the airframe, and varying the center of gravity on the natural fre- 
quency variation is shown in Figure 4.6. The variation of the natural 
frequencies for aerodynamic conditions 1 and 2 is independent of the dither 
frequency in the range from 10 to 40 radians per second. The variation for 
condition 1 is approximately 0.034 radians per second and for condition 2 
is approximately 0.123. The reference natural frequency for condition 1 
was 1.355 radians per second and for condition 2 was 1.185. The variation 
of the natural frequency depends on the dither frequency for aerodynamic 
conditions 3 and 4. The average variation of the natural frequency for 
condition 3 for dither frequencies in the dither range is approximately 
0.068 radians per second. At dither frequencies of approximately 23 and 
24 radians per second there are no variations. The average variation of 
the natural frequency for condition 4 is approximately 0.06 radians per 
second with a minimum variation of 0.028 at a dither frequency of 40 radians 
per second. The reference natural frequency for condition 3 was 1.225 
radians per second and for condition ^ was 1.51. 

In summary, the variation of the natural frequency does not depend on 
the dither frequency for aerodynamic conditions land 2 but does depend on 
the dither frequency for conditions 3 and 4. The variations were small for 
all conditions with 10 per cent of the corresponding reference value being 


the maximum variation determined, The effect of the flexibility of the 
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airframe and of a varying center of gravity on the variations of the 

natural frequency is not clearly indicated. The magnitudes of the average 
| variations of the natural frequencies for the aerodynamic conditions varied 
| from large to small in the following order: condition 2, 3, 4, and l. 

The investigation of the variations of the damping ratios and natural 
frequencies was conducted with the amplitude of the closed loop response 
selected to result in an average damping ratio during the flight of 0.7. 
Since the defined optimum performance of the system was a damping ratio of 
0.7 with a minimum variation of the natural frequency, the optimum dither 
frequency, or frequencies, for each of the aerodynamic conditions was 
determined on the basis of minimizing the variation in damping ratio and 
natural frequency. The error criterion used to determine the minimization 
of the variation in damping ratio and natural frequency was VAS + Aw2 . 
Plots of Ad versus AWw are shown in Figures 4.7 through 4.16 for four 
aerodynamic conditions. The optimum dither frequencies selected for the 
aerodynamic conditions are shown in Table XI with the corresponding error 
criterion value, the variation in the damping ratio, the variation in the 
natural frequency, and the Reference signal value, Two dither frequencies 
with different error criterion values are listed for condition 4, The dither 
frequency of 40 has the lower criterion value, but the Amplitude versus Kg 
curve for that frequency has a very small slope at the Reference value which 
would require a very accurate setting of the Reference signal. The dither 
frequency of 36 has a larger slope and would be the preferred frequency if 
the system was actually mechanized. 

The Reference signal values listed in Table XI are equal to the closed 


loop response amplitudes required for an average damping ratio of 0.7 for 
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Aerodynamic Opt imum Error Variation | Variation Re ference 
Condition Dither Criterion i in Signal 
Frequency Value Natural Value 
or Frequency 
Frequencies (kad/Sec,) 
(Rad/Sec. ) 
Condition 1 10 0.56 
40 0.42 
Condition 2 05 
0.73 
0.7 


0.375 0.027 0,993 
0,405 | 0.028 | 0.975 


72 











| 
| 





the ۲ The is e xréct 1f tri رت‎ hálr malt 816 ۶ 
the dither input siensl end the gains of the 于 RS and rectifier in the 
adaptive loop are equal to unity, [f the Jails are nct unity, the 
Reference value listed must be corrected by being multiplied by the pro- 


duct of the half cycle average of the filter, and the gain of the rectifier. 


4.4 Application of the Results to a General Case, 


The object of this section is to formulate a procedure from the empiri- 
cal results discussed in Section 4.3 and from the analytical verification 
contained in Appendix C to facilitate the rapid selection of the optimum 
dither frequency. The selection will be based on minimizing the damping 
ratio variation. Once the dither frequency or area of frequencies is 
selected, the procedures of Section 4,2 can be used for the frequency to 
determine the Reference value and the magnitude of the damping ratio varia- 
tion for the average damping ratio desired for the flight. 

The dither range is established by selecting the lower frequency limit 
five or more multiples greater than the maximum short period frequency to 
eliminate coupling. The maximum frequency of the short period mode can be 
approximated by the square root of the square of the real and imaginary 
parts of the zeros of the Inverse Model. The upper limit of the range is 
fixed by the actuator requirements at 40 radians per second, 

The assumption is made that the controi system will be investigated by 
a point study of the flight trajectory, From the root loci shown in Figure 
4.11, it is noted that the open loop poles associated with the bending modes 
remain relatively stationary during a flight but that the pen loop zeros 
associated with bending modes move during the flight from smaller to larger 
values, If slosh modes are temporarily ignored, che effect that the varying 
values of the bending zeros have on the shere perind mode is to vary the 
shape of the short period characteristic, The bending zer^s must be in the 
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ES ہاج رر‎ the ۶ 3 1 | n eppretiable effect 
bn tha د‎ ihe flisht hase with tro hand : w.tb the smallest 
imaginary pact, which is the earliest flight csi» —onsidered after launch, 


results in a short period charactertisti: ኮጊ...ከ hes 2 sarger adapative gain at 
a particular damping ratio than any of the other cases, The open loop poles 
and zeros associated with the short period appr.»'mitınn alsz vary during 
the flight, but it is assumed that the variations are negligible compared 
to the variations of the bending zeros. The values of adaptive gain at a 
particular damping ratio decreases as the imaginary part of the bending 
zero increases. The adaptive gains for a particular damping ratio are, 
therefore, inversely proportional to the time in flight. The dither system 
maintains the amplitude of the closed loop response constant during a flight. 
A constant amplitude results when the product of the adaptive gain and the 
dynamic gain remains constant. The dynamic gain is the gain associated with 
the vehicle dynamic characteristics, Tn order to have 2 minimum variation in 
the damping ratios, the product of the dynamic gain and the adaptive gain 
corresponding to the desired damping ratio must be approximately constant 
during the flight. Since the adaptive gain varies inversely with the time 
in flight, the dynamic gain must vary directly with the time in flight for 
the product to be a constant. If the dither frequency is selected close 
to a bending zero of the initial flight case, the dynamic gains will vary 
directly with the time in flight. The bending zero must have a value with 
a small real part and an imaginary part in the dither range. 

Tf the value of the bending zero has an imaginary part in the dither 
range but a large real part, the selection cf the dither frequency for 
optimum damping ratio variation is independent of the position of the zero. 


The value of the dither frequency does not appreciably affect the dynamic 


د : 








٧ |] "۳۱ ۱۲۳۷:۱۲۵۶ 1 ۱۷۱۵5 ۱ dompinp -! SOT 72888023. During 
the flight tbe bending modes do result in ۳/۵۵5۵ LT ihe adaptive gains 
for the same damping ratio, Since the val. ተ . Gicber frequency can 
not result in dynamic gains varying inversely tc the short period gains, 
large variations in the damping ratio may result. The requirement of deter- 
mining when the real part of the value of the bending zero is large or small 


jue of the dither frequency 


3 


can be eliminated by always selecting the v 


t 


close to the imaginary part of the bending zero of the initial flight case. 
If bending zeros do not exist in the dither range, the optimum damp- 
ing ratio variation is independent of the selection of the dither frequency. 
The value of the dither frequency can not vary tae dynamic gain to result 
in adaptive gains which correspond to the values required for constant short 
period damping ratios, The variations in the adaptive gains of the short 
period modes for the same damping ratio increase as the number of bending 
modes increase; and therefore, the magnitude of the damping ratio variation 
would be expected to increase as the number of bending modes increase. 
Slosh modes may have an important effect on the dynamic performance 
of the system for some airframes, Slosh modes would have a negligible 
effect during the initial and final phases ct the flight tecause the fuel 
tanks are then approximately full or empty. The slosh modes would affect 
the middle phase of the flight. [If the slosh mode frequencies are close 
to the short period frequencies, the shape of the short period characteristic 
could be appreciably varied for the middle flight cases, The shape of the 
characteristic would not be affected for the initial and final flight cases. 
The largest adaptive gain for a particular damping ratio would be for the 
middle flight case which is most affected by the slosh modes, The dynamic 


gain would have to be reduced by the selection of a dither frequency for 
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۲۱۳1 dinah ۱ I ۳ pes yst : 501۳ ۲ نه‎ large adaptive 

zair if a mall damping fário ۲۱/۱8 کا‎ ,5 "tai-ed. A dither 
frequency af b value zwual ta ihe zihar. i a modings zero of the 
middle case wıll reduce the dysamıc gain if the real part of the zero is 
small. The dither frequency selected, however should be between the values 
of the imaginary parts of the bending aero of the initial case and the bend- 


۹ 


ing zero of the middle case to result in a small demping ratio variation 
throughout the entire flight. The exact frequency in the area would have to 
be determined by the procedure in Section 4.2. 

In summary, the selection of the optimum dither frequency ís as follows: 

1。 If no bending zeros exist in the dither range, the value of the 
dither frequency is not restricted, 

2. If there are no important slosh modes in the system but if 
there are bending zeros in the dither range, the optimum frequency is close 
to the imaginary part of a zero associated with the initial flight case, 

3. If there are important slosh modes in the system and there are 
bending zeros in the dither range, the optimum frequency area is between the 
value of the imaginary part of an initial flight case zero and the imaginary 
part of the zero of the middle case which is most affected by the slosh modes. 

The dither frequency must not assume a value that is close to the value 
of the imaginary part of a pole which is adjacent to the imaginary axis be- 
cause the adaptive gain change per unit amplitude of the closed loop re- 
sponse is large. A large gain to amplitude ratio reguires a very accurate 


setting of the Reference value, 
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E HA F ፈ E R ጋ 
CONCLUSTONS 


The conclusions from the results of the investigation are as follows: 

l. The dither adaptive loop employs a single parameter to control 
basically a single dynamic characteristic of the system. The parameter of 
the loop is the closed loop response amplitude and the characteristic con- 
trolled is the damping ratio. 

2, The components of the loop are simple to mechanize, but a prior 
knowledge of the airframe dynamics are required to select the dither fre- 
quency, Reference value, and the shape of the gain curve for the adaptive 
gain pot. 

3. The adaptive loop will adjust the adaptive gain to partially 
correct for the varying dynamic characteristics of a flexible airframe if 
the dither frequency has no coupling with a bending mode, The accuracy 
of the correction in maintaining a constant damping ratio is the inherent 
accuracy of the Amplitude Closed Loop Response Criterion. 

4. The adaptive loop makes satisfactory corrections of the adap- 
tive gain for disturbances. A further investigation into the gain curve 
shape for the adaptive gain pot and into the optimum locp gain is required 
to verify that an actuator signal will not cause an unsetisfactory variation 
in the adaptive gain. 

5. The Amplitude Closed Loop Response Criterion does not result 
in a constant damping ratio. The variation of the damping ratios increases 
as the desired average damping ratio decreases, as the airframe flexibility 


increases, and as the variation of the center of gravity increases, If 
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frequencies in the dither range greatly affect the amplitude of the vehicle 
dynamic response, the variations of the damping ratios can be appreciably 
reduced by the proper selection of the dither frequency. 

6. When the damping ratio variation varies with the dither fre- 
quency, the minimum natural frequency variation occurs at approximately the 


same dither frequency as the minimum damping ratio variation. 
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APPENDIX Á 
The aerodynamic equations for the missile were chtained from the 
Autonetics Division of North American Aviation. The equations are two 
degree of freedom short period mode approximations with three bending 


modes and one fuel siosh mode included. 


"۰-۰ عم‎ 6 1C,8 oh + Eid) + Fi alls + Elda + me + 
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wher 
ٹوچ‎ angle oí attack in degrees 


= 7ک‎ pitch angle in degrees 


dı er first bending mode in feet 

da = second bending mode in feet 

d3 - third bending mode in feet 

Se - pitch angle as measured by position gyro 
UNT — first ignition fuel slosh mode in feet 


/۱ )رخ‎ 一 first liquid oxygen slosh mode in feet 
EH — first liquid hydrogen slosh mode in feet 


ST == thrust control deflection in degrees. 


t-20 t-75 t-101.9 t-200 
Ay - „064 - „0489 - „0336 - „0236 
B, 1 1 1 1 
of - „158 - „0226 - „0103 „0074 
0 0 „00123 „000229 .00015 
E, 0 „0337 „031 - „061 
F, 0 - 000886 - 000167 -,00019 
E 0 - „0461 - „0154 - „185 
H, 0 001108 ‚000214 . 00014 
7 0 „0940 „0346 - „047 
3 „0567 „0289 0 0 
Ki „058 „0298 „0422 - „039 
P. .00497 „00255 .00362 - „0027 
Mi „222 „0565 „0378 „02 
E 0 „0217 „0159 1.65 
Ba 0 -,00156 -, 000519 -,038 
جک‎ 0 - „00204 -.000813 -,055 


o? 








ሙሬ!) 2-75 ESO t- 200 
ps 0 0515 - 0386 -31,92 
E: 0 00168 . 000664 „042 
Fa 0 „071 „0552 -12.73 
| Ga 0 „00213 -, 000792 -,024 
| Ha 0 - „163 - „127 -28.33 
| en „00316 „248 0 0 
| is -, 0048 -,00716 .00193 -3.8 
Ko - 000919 -, 00205 -, 00482 .012 
| La .0485 .0814 E 2,28 
A3 0 401 „213 32 
| B 3 0 - „0294 23ء ۔‎ - „012 
af -.163 -.205 -.175 -.59 
Dz -66.5 -67,3 -66.9 -767.46 
E 0 „0348 „00981 „023 
Ez 0 1.39 „821 16,34 
Ga 0 - „0482 - „0125 -,0013 
Ha 0 -.3.24 -1.83 6.0 
i 137 - „0414 0 0 
m „131 „298 „634 3.08 
K - „0001 - „0182 - 9 - „128 
7 -1. : -1.79 -1.86 
Ay 0 - „288 - „155 - 9 
By 0 „0243 .00591 .009 
Cy 0 „0348 .00981 „023 
Dy 0 . 708 „384 13.82 
Ey - .43 - „463 - „439 - „96 
Fj 2 -463 -463 -2169.63 
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Lg 


Aa 
A3 


The magnitudes of the remaining coefficients were 


included in the transfer functions. 
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0 
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„0544 
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„199 
- „00704 
- „0125 
M 
„0.24 
„584 
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‚015 
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not available but are 








ቆዕ e 
' ١ LA Y exuations @ and X will be shown. Information 


ef 


ES‏ دو e‏ وه وه 
Ay j ZA x , and Ag‏ و (MOOV fhe jistification of the dy , de , dz‏ 9 


4 


equations 1s contained in Reference 10. 


From Newton's second law of motion. 


EE, - Ale “መ 
گور ری یں‎ 
A 


ot LE 
oy ELF Dz +.(oxh) 
E 1 Zl ማፕ ጭጭ. 
£M= áa í 
EN = har + £-(WXh 


where a, 2 E ; 2 Es, are the forces in the x, y, z directions of a 


right hand system of Cartesian axes fixed in space; where £L, £ M, £N 
are the morents about the x, y, and z axes; where hæ ۱ = , and nos are 
the moments of momentum about x, , y, , and z, axes which are fixed to 
the airframe; where xh), Gh) , and AX ) arise from 
the angular velocity W of the Ky والا و‎ and z, axes with respect to 
MA ول‎ z axes; where ax, 2 y = سوه‎ , are the accelerations in the 
RE and 2 ገ... and where m is the mass of the airframe. 
The absolute acceleration of a body with respect to x, y, z axes is 
A ص‎ 
aale = IE TWXV 
W 


where V is the instantaneous linear velocity and is the angular velocity, 


Ly - (۷ ٧ 
W XV =(QW-RvV)e + (RU- PW Ji + (PY-Qu) ሳሎ 
where F, Q, and R are the angular velocities along the x, y, and z axes 


respectively and U, V, and W are the linear velocities along the x, y, and 


z axes. The components of the acceleration are then 
Qe ae = Ur QW-RV 


5ھ 








en - V + RU- PW 
24, 2 W tPV-QU 
Substituting the components of acceleration into Equatíons 1 and consider- 
ing the mass as a constant the equations become 
EF, > (Ur qQw-RV) 
2 £Fy = و لاگ‎ - Pw) 
2 Fz = am (WEPV- Qu) 


The components of the moment of momentum are calculated from summing the 


moments of the velocity vectors about each axis and multiplying by the mass dm, 
ha = Play - GL, - Rigg 

ካሃ TF 9 ማያ 

ha : له‎ Ay DA 

where Lo = (yt 32) ore and Tøy = 4ራ "Lm 


XY and XZ are assumed planes of symmetry, and therefore 


Equations 4 become: 


Ec = Ue = و اا‎ O. 
ከፉ = PI g% 
(5) : 
ህይል. 
ha = Rl ag 
Les = Ba 


£y: 9T, 


= እ = Als . 


(6) 


Since 
ui x) s (ha Qo ho R)e Oe R- hs Pg Pho p) 


Equations 2 with Equations 5 and 6 substituted become 


— 
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yl‏ و و97 سر ھا 
EM= Olay t PR Izy -Tyg-)‏ 


(7) ۰ 
l N= Ry y+ POT - Top) ٠ 


The forces indicated by the left terms of Equations 3 represent the 
summation of aerodynamic, thrust, and gravity forces. The gravity forces 


are shown below and their derivation can be found in Reference 9. 


ያራ E m pain Do) cow Dto Y + Img we Oo cere Po) 


coa- O sinr F- (mng cos Go cow Yo) سه‎ 


۴, = 02 p ہہ و‎ IV d) 
PE ^ 
771 (ሥታ رټسم‎ Oo an Po) (tow ሦ Ce D T av Poir- O oes i) 
(mg wen tr Jl ir © tn J) 


Egy Amin ltr at o D a Join PA 
ma al? Don Dr ao Vom 9 we P- re ford )+ 
(ong Cae Cr Cow J.) (or O Ceo 2 : 


where m is the mass cí the body; g is the acceleration of gravity; O, p, 
and y/ are the Eulerian Angles for a ZYX system and 6 ره‎ Do, and Wo are 
the steady flight values, 
/ 
Defining Era ۱ 2 y و‎ ۰ £ FA as the sum of the 
aerodynamic and thrust forces, Equations 3 can be written as 
=, . 
ደ/ሠ = me (Ut QW- RV) - Fog 
8 7 , 
> 2 Fog on(V +RU- PW) -ፖ 


#8 = m(wtPV-Qu) - Fag 


17 
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Each of the i :stintaneous velocity components can be written as the 


sum of a steady state value and a change in velocity. 


U: ولا‎ + 1e- 
v z Vo ተ” 
W: Woru 
ዌን. tp 
Q = Qo rg 


h* Bera 


ሠ 2 ወ ጋ. 





In order to linearize the force and moment equations it is assumed 

that 

1) the velocity changes are small and that the products and squares 
of the velocity changes are negligible compared to the velocity changes 
themselves 

2) the disturoance angles are small and the sine of the angles 
can be set equal to the angles and the cosine of the angles set equal to 
one. The product of the angles is zero, 


Equations 7 and 9 can now be written as: 


EFE solae QoWo +Wog + Qour-RoVo=Ronr- Vos + 
ger Bo)- en bn at PY + lens com JO] 
(10) 
EF) E رل‎ 0 mp 
2ص و‎ B LP] 


£F flot Rr ولاه‎ DA 
سار"‎ N Did y) 


(11) £L C Epp FS Re T Qu act Rog)l T33 -Igy 
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2 M= ያ ] /ሃ t (2 Kat E Rage) LT “12 7 


EN = رز رآ‎ FQ Ru ره کال ام‎ Teg. 
The equations which relate angular velocities about the x, y, z axes 
to the Eulerian Angles ©, 7 V are 
Pz 0 -y un” E 
Q= we Pr Yan Pin © 
R= ¥ " می‎ P tor 9 - Oo d. 


When linearized the equaticns become 


p= 
(12) 9:8 
- ې‎ 


The steady flight condition is defined as all velocity components equal 


to zero except Up and Wo. Using this definition and Equations 12, Equa- 


tions 10 and 11 become 

EF) = an Lat + Wag ور‎ tne Qo + g (cee Bo) ره‎ 

mart - "Wo. - ge 9)V- o 92g]‏ < بوک 
E <arLir=Uo gtg (et (69 4 cee Gol‏ 2 


(13) 


)14( 


۱1 


27 PL, 
íM=ð Ij4 
Nee موی وتا‎ 


The aerodynamic forces and moments can be represented by a Taylor 


4 


series expansion. Because of the planes of symmetry and the steady flight 
Y ሥ : 
conditions ሩ ቸራ : کت‎ , and M are functions of only q, u, and w of the 


velocity components and Fy , L, and N are functions of only p, v, and r, 


39 








chis |S proved fn ۱ ۰:۵6 @, The other variables affecting the forces 

and mmwıts are the angie of thrust deflection, the airframe flexibility, and 
the changing mass of the body. Since the changes are assumed small the 
second derivative and above for the velocity and bending terms are consider- 


ed negligible and the first derivative and above is considered negligible for 


the change in mass, The resulting equations are of the form 


oy Bae ra ria oda مرا رر‎ 
+ re) Ws GE abs. T tee rt ay 
“ረ. + (Ez, گے(‎ + HIR 2 us (#2) da 


3 (Brox 十 Gr) A 10x د‎ Gk, A LH 


where d, , 1, , and ds are the first, second, and third bending modes; 


PA 


or is the angle of thrust deflection; and AL ; 人 LOX , and ALH are 
the changes due to the change in the mass of the fuels, 
The thrust axis is assumed to pass through the center of gravity and 
hence does not contribute to the moment equations, The thrust force equa- 


tions are 


£ Fa BIR 
= E, = -Tain é 


T 
where 3 is the angle between the X axis and thrust line. 


۱ ۰ AT 
and 4E z- 1 


Therefore 


ca cov E. dE o‏ ديم 


(16) سب‎ 
Fy = - E E A 


= IN a (Vs. 
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To determine the equations of steady flight substitute Equations 16 
and 15 into Equations 13 and 14, substitute the values of the steady flight 
condition, and set the velocity changes equal to zero. The resulting equa- 


tions are: 


Fra + To oe Ezme grains Do = 


Fago 2 2 
(17) i 
Fare IR tg مسن‎ = e 
lO 
Mo =O 
۸۷ 2 ۵ , 


When stability axes are used as the reference system and quasi-steady 
flow conditions are assumed, all terms containing W, disappear and all 
aerodynamic partial derivatives with respect to rates of change of velocit- 
ies and with respect to rate of change of angle of thrust deflection are 
eliminated, For the Two Degree of Freedom Short Period Mode Approximation 
هم‎ is set equal to zero and only the equations ur and ZM are 
involved, Substituting Equations 16 and 15 into Equations 13 and 14 and 
subtracting Equations 17 gives the following equations of motion for the 
disturbed body: 

ሪሪ”ሥ፦ Ve gr Pg Gi Go )O = ይሪ 5. + Fare tE #7 Sr + Ff, ተ 


ELL, +ያ رك‎ + da + Fáa ds t 
ይፊ 2: + FAT Ar + Farox OLoxt Farn ALH 


(18) = 
A Les + Ma, Ka 

LAWA t Es Te Me: 4 + Mag AL + Maron ALOX Ts 
MALH ALH 


JE T M, 
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the an: ۲ ۱ 3W ہہ‎ rs de fined as the angle between the wing 


chord line S- tha r. lt ۷/۰۰ wind, 
7 ላ/ fie ፡ خي‎ ^w 
A -ሙ سے وو گے‎ 
Ad “ላጪሟሚ። [/ Ue y rte Uo SNCE MEÐ 
orn ۳ 
(19) nal ےی‎ € Ubpot. 


Equations 18v::h “quations 12 and 19 substituted become: 


à Fara t nie) o ra) y (E y, 


60 +) Ê) da + Ls + (Bahr 
(20) (ÆÐ Ár (E829) A Lx + (EE) A LH + (FE) Sr 


@ =(MurUs oe + (M) 2 +(Md)L + (Md Jdi + ML) La 
FM ) ሪሬ ተ (miei) s + (Mts) s+ (Mar) AT 
+ (Marox )ALOX +(Mark) OLH , 


The coefficients for any instant are constants. Equations 20 agree in 
> n 
form with the given equations for ol and ©. 
The Transfer Functions which were obtained from the Aerodynamic Equa- 


tions are as Follows: 


1. Aerodynamic Condition 1: 


NN م9‎ - 6 ee t. 064) 
v a ¿lot 064) 


Was 75 Op = NH K^ (71.068) 
ár A (a- ME) +02) 


Ez/01.9 مه‎ - Job (st.0374) 


AT ورا‎ ٧ ጋታ 7፡2 


...ے38 +ے) 2.28 
ST‏ 


(¿0074 (6-4 42 (24 4.32.) 





۱ 
دل‎ 
0 
0 
(D 
"S 
M 
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Ee.‏ مور 
p^ ST‏ 
1T. 0 -6,06‏ 27[ 032, رت ol‏ 





Sr Au. lar LE 


pue 200 ANS Co 777, 22 
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3. Aerodynamic ۵۵۲ 41610۳8 3: 





dos” GULo-T. o4 6 t 917‏ ۶ ص)( 8 قص تصا اه ez‏ :| دعر 
A‏ 742 1 


ر2 2 D — ተ. MA‏ الع د 














uox 


Cp 2 i 
p ALO ረር ee = 
PT ዕ ‹‹‹“፲/ዕ-ሪጋሃ/ኃታሪ «/ሬታ. 2272 ty angor 7یت‎ 74277311) 
4. Aerodynamic Condition 4; 





EP = 10.78(ar:06‏ وو دغر 
F ST E A, hof Tas): 917 (21.027 14 27‏ 








Hol 9 Be ~ کش نک‎ 1 Io +972 ھت‎ 2) 
ár  (a-el)G- E 73)(et. ÆÐ ar ,#5/ 


Car › ዕ3ሬ ! >.// ኃሟያ(ፊኔ t.06 6.18069 r.032. ty 3.33) 
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ANALYIECAL VERIZTCATIORN OF Wee ባዊ SLA ه٢ و اسا ده‎ 


The results of the variation of damping ratios with dither frequency 
obtained by empiri al methods in Chapter 4 will be analyzed ty analytical 
methods in this Appendix. 

For the analysis the root loci for the flight cases of a particular 
aerodynamic cord tion will be assumed to have the same shape and system 
gain values, The transfer functions of the vehicle dynamics for a parti- 
cular aerodynamic condition do not vary to the extent that the use of the 
assumptior to determine general trends would be prevented. A particular 
syster gain, K, on the roct loci for the three flight cases of a condition 
would determine the same demfing rátio,á . The system gain is the produc: 


9 


of the dynamic gain, Mg , and the adaptive gain, K$ . The dynamic gain 


pt 


pie 
ፍን 


the gain associated with the vehicle transfer function. Using condition 
3 as an example, the dynamic gains for the flight cases are approximately 


equal to the following values: 


Case t - 20 Ms = 16.05 
Case t = 75 | 
Case t 101.9 Ms 21.6 


Since a particular damping ratio has the same value of K for the three 


cases the values of Ks for the cases are: 


RS ze ب‎ ۳ — 
ES 415 1፡28 


K § ሯ= /6ፃ - 216 
and Ket & 20 70ا یج‎ > Kg t - 101.9 for a particular ፊ . 


plot mi S versus ys is of the following form. 
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Ks 


The trends shown in the á versus Kf plot are valid for any of the four 
aerodynamic conditions, 
The closed loop transfer function for the system is 
هه‎ a Sr 
سے ےنت‎ ni T 
"۳000۸ ۶ھ‎ It K> ^ھ)‎ + 3.6 3.6) E 
as 
w^en a is the vehicle dynamic transfer function. When the dither freq- 
uency is large compared to the poles and zeros of the vehicle dynamic trans- 
PA 
fer function, the transfer function is approximated by QE = m 
ST መራ 
This approximation is justified for the dither range, 10 to 40 radians per 


second, for conditions 1 and 2. The closed loop transfer function becomes 
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Since the dither system maintains the amplitude cf the closed loop response 
constant, 


Ks Ms = KaM - KsaMs2 ات‎ 


and for Aerodynamic vonuition 1 
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and for Aerodynamic Condition 2 


Ks, 8 onn. [Mc 9 
b+ Ks, S18 | ተ/ዕፅ ሬ፡ፖ | ۲5 2 


The relationships are independent of the dither frequency, and constant 
damping ratio variation would be expected in the dither range. This agrees 
with the empirical results for conditions 1 and 2 as are shown in Figure 
4.4. ፲፻ Kg is assigned the value 1 for both condition 1 and 2, then for 
۲٣۰۰۰1۴٦165 1 

ሾና | = 1,67 K $3 = 07520 
and for condition 2 

> 1, = 0.72 

KS 1 B. K $3 U 
For both conditions, Kg 1? KS 2 2 K S 3 which ís the trend necessary 
for small damping rati. variation as determined from the d versus Kg Dior, 

For condition 3 the dither frequency can be close to poles and zeros 
of the vehicle dynamic transfer function. Tf the dither frequency ís close 
| or . 

to a pole of the transfer function, then — is very large and 
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and 
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The amplitude of the closed loop response is independent of the adaptive 
gain, and therefore, Kg could assume a infinite number of values with 
an infinite number of damping ratio variations, 


If the dither frequency is close to a zero of Case t = 20, then 
E. WO ረ E A 
OTL Jo en 757 STL= lol. 7 
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; amplitude 3f rhe closed loop response, Ks Loos) 


= 1 





) Ket መ 25° ፈፈ t 101.9 which is the trend for a small damping ratio 
/ 
variation. If the dither frequency is close to a zero of Case t = 101.9, 
then 
መ. y Or 
مه رل‎ STL=75 STA ze)? 


In order to have the same amplitude of the closed loop response Ks t 207 

Ks | 2 و‎ t = 101,9 which is tke trend which would resoit in a 
large damping ratio variation, If the dither frequency is close to a zero 
of Case t = 75, the magnitude of the damping ratio variation would vary 
between the magnitudes when the frequency is close to a zero of Case t - 20 
and when the frequency is close to a zero of Case t = 101.9. 

The damping ratio variation curve for condition 3 should, therefore, 
have an undetermined damping ratio variation at the pole at approximately 
40 radians; have a small value for the variation at the zero of case t = 20 
at approximately 24 radians; and have a large value for the variation at 
the zero of case t = 101.9 at approximately 27 radians. This agrees with 
the condition 3 curve shown in Figure 4,4, 

In order to determine the dither frequency that the minimum variation 
in damping ratio occurs, the ideal expression for no variation which is 


of the form 
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and for a damping ratio of 0.7 


VF 


NS t = 20 = 0.86 


K ہم‎ 二 0.26 


Qn 


TOT 58 


If 


ኪሩ t 


is evaluated for various dither frequencies close to the frequency of 
the zero for case t = 20, and the smallest amplitude variation between the 
three cases sould indicate the optimum frequency. A dither frequency of 
23 radians per second resulted in the smallest variation in the amplitude. 

For aerodynamic condition ^ there are no zeros or poles of the vehicle 
dynamic transfer function which are close to the imaginary axis in the 
dither range.  Varying the dither frequency does not result in wide varia- 
tions in the values of the vehicle transfer functions, and the closed loop 


transfer function can be approximated by 


© م 


an Ww WE. 
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A constant damping ratio variation would be expected and agrees with the 
curve for condition 4 in Figure 4,4. Since the values of the vehicle 
transfer functions can not be varied appreciably, the variation in the 
adaptive gain required for a constant damping ratio can not be realized 


and a large variation in the damping ratio results. 
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